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Abstract 

We are interested in electrons kinetics in a stellar atmosphere to validate 
on : or invalidate the usually accepted hypothesis of thermalisation of electrons. 

For this purpose, we calculate the velocity distribution function of electrons 
by solving the kinetic equation of these particles together with the equa- 
tions of radiative transfer and statistical equilibrium. We note that this 
distribution can deviate strongly from a Maxwell-Boltzmann distribution if 
non-LTE effects are important. Some results and astrophysical consequences 
"■■-^ ■ are examined. 

O ' 1 Introduction 

c/5 Our work is concerned with the kinetics of electrons in a stellar atmosphere, mod- 
^ ' elled as a parallel-plane slab irradiated on a face. Our models of atmospheres start 
in the deep layers of stars, where the radiative field can be described in the diffu- 
k> sion approximation, and end with the layers of minimal temperature, before the 
. chromospheric raise whose effects are ignored. The free electrons are characterized 
' by their velocity distribution function : the electron distribution function (edf), 
which is calculated with the other thermodynamical quantities of the atmosphere. 
Our main objective is to understand the mechanism leading to the thermalisation 
of electrons, where the edf tends toward the Maxwell-Boltzmann distribution. It is 
accepted, in stellar atmospheres theory, that the thermalisation of electrons is ef- 
fective as long as elastic collisions dominate inelastic interactions of electrons with 
the plasma, a rather well verified hypothesis for electrons having energies greatly 
below the first excitation energies of atoms and ions composing the atmosphere. 
This hypothesis is not necessarily correct for faster electrons. Our work follows the 
line drawn by some plasma physicists at the beginning of the 70s (Peyraud 1968, 
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1970; Peyraud et al. 1969; Oxenius 1970, 1974; Shoub 1977). Their work demon- 
strated the important role played by inelastic (collisional or radiative) processes 
in the equilibrium reached by electrons, which can deviate considerably from the 
maxwellian equilibrium at high energies. We present below a stellar atmosphere 
model which is not in local thermodynamical equilibrium (LTE), confirming the 
results anticipated in the 70s on the basis of mainly theoretical developments. 

2 The model 

This problem consists in solving the equations generally used to model a non-LTE 
stellar atmosphere (equation of radiative transfer, equations of statistical equilib- 
rium, pressure equation, equation of energy, conservation of charge) , coupled with 
the kinetic equation of electrons. This non-linear system is difficult to solve nu- 
merically because it contains two coupled kinetic equations : one for photons, the 
other for electrons. Therefore we have used the simplest model of non-LTE atmo- 
sphere : homogeneous (constant density of heavy particles n ), isotherm (constant 
temperature T), and composed with hydrogen atoms with only two energy levels. 
The deviation from LTE is then due to the escape of photons by the free surface. 
On the other hand we have included in our model the main collision processes 
existing in a stellar atmosphere (elastic collisions, collisional or radiative inelastic 
interactions). The elastic collision term of the kinetic equation of electrons is writ- 
ten in a BGK model with a velocity dependent collision frequency. This model 
accuratly fits the main properties of the usual Landau term (Fokker-Planck) . To 
solve the equation of radiative transfer, we used the codes of the transfer group in 
CRAL (Rutily 1992). 

3 Results 

In our model, we choosed n = 1.2xl0 23 m 3 and T = 6500 K, which are typical 
values in the solar photosphere. The plasma is optically thick at all frequencies 
(optical thickness greater than 100), leading to a high geometrical thickness Z 
since there is no temperature or heavy particles density gradient. Finally the 
atmosphere is irradiated on its internal boundary layer by a Planck radiation of 
temperature T. 

The figure la is a classical diagram showing the superficial regions where the 
non-LTE effects are important. Figure lb shows that the edf is not a Maxwcll- 
Boltzmann distribution in the non-LTE region of the atmosphere (see u/wr = 2), 
the deviation from a maxwellian distribution being important very close to the 
surface (z/Z < 10~ 15 , corresponding to an optical depth t < 30 in the Lya 
spectral line). In figure lc, we drawn the superficial edf at z — as a function 
of the electronic velocity. The edf tail of fast electrons is strongly depleted when 
electron energies are greater than the minimum excitation energy of the hydrogen 
atom (E = 10.2 cV, E/E R = 0.75, v/vb, — 0.866). The edf tail shows successive 
platforms centered on v/vr — 1.25 and v/vr = 1.65. 

These features were already described by the authors at the origin of this work, 
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referenced at the beginning of this article. The mechanism responsible for this 
effect is very well explained in Oxenius's monograph (1986), where the author 
outlines an interesting feedback effect tending to amplify the deviation of the edf 
from a maxwellian distribution. This mechanism starts when elastic and inelastic 
collision frequencies become comparable at high electronic velocities, which is the 
case for a weak ionization degree. 
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Fig. 1. a) Deviations of the populations from their LTE values as a function of the 
reduced geometrical depth z/Z. Coefficients bi(z) = m(z)/nf (z), where nf \z) are the 
Saha densities of the hydrogen atom in energy states i = 1,2, are used to characterize 
non-LTE regions (bi(z) ^ 1). b) Deviation of the edf to the maxwellian distribution 
f^{v) as a function of the reduced geometrical depth z/Z. Both curves are drawn for a 
given velocity v/vn — 1, 2, where vb. is the velocity corresponding to the Rydberg energy 
En — 13.6 eV. c) Deviation of the edf to the maxwellian distribution as a function of the 
electronic velocity v, at the surface of the atmosphere z = 0. 
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4 Astrophysical applications and conclusions 

Astrophysical consequences of this work are numerous. In general the deviation 
of the cdf from a maxwellian distribution has a direct effect on all thermodynamic 
quantities involving the edf, e.g. collisional transition rates or spectral lines pro- 
files. It has an indirect effect on all other characteristics of the atmosphere, because 
of the coupling of all equations. Transition rates are used to solve the equations of 
statistical equilibrium, which lead to the populations and ionization degree of the 
atmosphere. Inversion techniques of spectral lines observed by spectroscopy are 
also sensitive to the edf shape, so that temperatures or densities calculated with 
these techniques are a priori affected by deviation of the edf from a maxwellian 
distribution (cf. Shoub 1983, Owocki et al. 1983, Salzmann et al. 1995). Finally 
the non thermodynamical equilibrium of electrons may be at the origin of physi- 
cal processes which are still not very well understood at present, for example the 
heating of the sun corona (cf. Scudder 1992,1994). 

The results presented in this article are based on a very simple atmosphere model, 
which guarants a numerically stable solution. In non-LTE regions close to the 
surface, the edf shows important deviations from the maxwellian distribution in 
the fast electrons tail. Our model confirm, by means of numerical codes accurate 
enough to handle this complex problem, most of the physical ideas advanced thirty 
years ago. Our main contribution consists in the construction of a selfconsistent 
model of a stellar atmosphere with a priori non thermalized electrons. Also we 
have used very accurate radiative transfer codes. It remains to make this model 
more realistic for comparison with observations. 
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